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Nomenclature 


Throughout this paper Ti = c= 1. 

m = 938 MeV is the nucleon mass, which is used for both protons and neutrons. 
Eh is the total energy of the projectile (particle 1) in the lab frame. 

Tu is the kinetic energy of the projectile (particle 1) in the lab frame. 

Pn is the momentum of the projectile (particle 1) in the lab frame. 

10 3 GeV = TeV = 10 12 eV (Tera-eV). 

10 6 GeV = PeV = 10 15 eV (Peta-eV). 

10 9 GeV = EeV = 10 18 eV (Exa-eV). 

10 12 GeV = ZeV = 10 21 eV (Zetta-eV). 
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Abstract 


The total proton-proton and neutron-proton cross sections currently used in 
the transport codes HZETRN and BRYNTRN show some disagreement with 
experiment in the GeV and EeV energy ranges. The GeV range is near the 
region of maximum cosmic ray intensity. It is useful to correct these cross 
sections, so that predictions of space radiation environments will be more ac- 
curate. Parameterizations of nucleon-nucleon total cross sections are developed 
which are accurate over the entire energy range of the cosmic ray spectrum. 


1 Introduction 


A current assignment is to improve the elastic neutron differential cross sections currently 
used in HZETRN [1] and BRYNTRN [2, 3]. These cross sections and the relevant theory 
are described in Chapter 4 of Reference [4], where it is shown that nucleon-nucleon total 
cross sections are fundamental inputs to the formulas for the neutron elastic differential 
cross sections. It was therefore decided to check how well the nucleon-nucleon total cross 
sections currently used in HZETRN compare with experiment. It will be seen below 
that although the global fits are quite good, there are some discrepancies in the region of 
maximum intensity in the cosmic ray spectrum. The aim of the present paper is to make 
improvements. 

Note that the cross sections discussed in this work have also been parameterized by 
Kalinovskii et al [5]. 

The total (elastic plus inelastic) cross sections for proton-proton (pp), neutron-neutron 
(nn) and neutron-proton ( np) scattering have been parametrized in reference [4], The pp 
cross section is 


a. 


HZETRN 


pp 




a, 


HZETRN 


PP 


m, 



|40 + 109 cos(0.199-y/rh) exp 


exp < 6.51 


exp 



— 0.451(7b - 25) 0 ' 258 
(for T u > 25 MeV) 
(for Tu <25 MeV) 


( 1 ) 

( 2 ) 


Note that this second equation is written with a typographical error in reference [4], with 
the minus sign in the wrong place. Note however that this typographical error is not the 
cause of the discrepancies noted in the present work. The actual formulas currently used 
in the transport code HZETRN [1] do not have the typographical error and are written 
correctly, but nevertheless show some disagreement with experiment as explained below. 

The total cross section for neutron-neutron scattering is not available and therefore 
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the proton-proton cross section is used, namely 


Vnn(Tll) = <Tpp{Tu) (3) 

This is accurate due to isospin symmetry. The neutron-proton cross section is 

<f™(T u ) = 38 + 12500 exp [-1.187(Tu - O.l) 0 - 35 ] (for T u > 0.1 MeV) (4) 
a HZETR. N (Tu) = 2 6000exp -(T^)° 3 (for Tu <0.1 MeV) (5) 

Again, this second equation is written with a typographical error in reference [4], with 
the minus sign in the wrong place. In the above equations, Tu is the projectile kinetic 
energy as measured in the lab (target) frame in units of MeV. The cross sections are in 
units of mb. 

Many parametrizations found in the literature, as well as experimental data [6], are in 
units of GeV and therefore, when comparing to data and other work it is more convenient 
to use GeV. The formulas above take on a different form when expressed in GeV. In order 
to do this note the following. Consider equation (1). The terms 5 /Tu and the arguments 
of cos and exp are obviously dimensionless. This means for example, that the 5 should 
really read 5 MeV, with the same units as Tu. Thus, in GeV this must be written as 
5000/Th. Also the number 0.199 must have units of energy -1 / 2 . The number 25 must 
have units of energy, and in GeV units it becomes 25,000 GeV. The number 0.451 must 
have units of energy -0 - 258 , so that it becomes 0.451/1000 0 - 258 GeV -0 - 258 . 

The formulas are conveniently re-written as 

a^ET RN (T n ) = (l + 1 40 + 109 cos (b^) exp [-c(T u - d) 0 ' 258 ] j 

(for T lt >25 MeV) (6) 

Vrr ETRN (Tu) = expje.51 exp^-^)^ j (for T u < 25 MeV) (7) 

<Jn p ZETRN ( T u) = 38 + 12500 exp [— /(Th — g) 0 ' 35 ] (for T n > 0.1 MeV) (8) 

<? ETRN (Tu) = 26000 exp ~ (for T u < 0.1 MeV) (9) 

where the constants a, b, c, d, e, /, g and h are given in Table 1 in the MeV unit system. 
It then becomes a straightforward matter to convert to GeV. 
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Table 1: Constants for Cross Section Formulas (6) - (9). 


a 

b 

c 

d 

5 MeV 

0.199 MeV" 1 / 2 

0.451 MeV” 0 ' 258 

25 MeV 

e 

f 

g 

h 

134 MeV 

1.187 MeV -0 ' 35 

0.1 MeV 

0.282 MeV 


2 Comparison of HZETRN formulas with experiment 

The HZETRN and BRYNTRN formulas for proton-proton total cross sections are shown 
as the dashed blue curves in figures 1 and 2 compared to experimental data, shown as 
the solid triangles, from the Particle Data Group [6]. Figure 1 shows that the apparent 
overall agreement between theory and experiment is quite good except in the EeV region, 
where the cosmic ray intensity is very low and therefore of negligible concern for space 
radiation. The GeV region is of most concern because it is near this region where the 
cosmic ray spectrum reaches a peak. Therefore, this region is examined in more detail 
in figure 2. The HZETRN cross section is found to oscillate which is not good behavior. 
In fact, this oscillation is present at most energies. Furthermore, the HZETRN cross 
section shows some disagreement with experiment in this energy range. Comparison with 
neutron-proton data is shown in figures 3 and 4. Again, the global fit looks quite good as 
shown in figure 3. However, when one looks more carefully at the important GeV region 
as shown in figure 4, one again notes some discrepancies. The aim of the present work is 
to correct these discrepancies. Their significance is that transport codes may give slightly 
incorrect results. The input cross sections should always be as accurate as possible. 

3 New parametrizations of nucleon-nucleon total cross 
sections 

The GeV region is often referred to as the intermediate energy region. The high energy 
region is up to 10 9 GeV. Recall that 10 3 GeV = TeV = 10 12 eV and 10 6 GeV = PeV = 
10 15 eV and 10 9 GeV = EeV = 10 18 eV and 10 12 GeV = ZeV = 10 21 eV. Thus, the region 
being referred to as the EeV region is of great interest from a pure physics point of view 
but is not of much interest for space radiation because the cosmic ray flux is very low in 
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this region. Nevertheless, for the sake of completeness, parametrizations in this region are 
also studied because data are available. The dominant process in the very high energy 
region is due to Pomeron exchange, such as flavorless gluons [7]. Parameterizations of 
this very high energy region has been completed by Collins and Martin [7] , Wong [8] and 
De et al [9] and many other authors. This work has been condensed into a formula that 
provides the best fit between intermediate and high energies by the Particle Data Group 
[6]. Their cross section parameterization for both pp and pn reactions is 

a PZX? (7n) = Z ab + Blog 2 (s/s 0 ) + Y 1 a \s 1 /s)^-Y 2 ab (s 1 /sr (10) 

with the square of the center-of-momentum energy given by 

s = 2 m(m + \Jp\i + m 2 ) ( 11 ) 

where m is the nucleon mass and 

pu = + 2 m) ( 12 ) 

In equation (10), the constants Z ab , B, Y( jh , Y/f 1 are in units of mb and s, s i, So are in units 
of GeV 2 . The values of the constants are listed in Table 2. 

Table 2: Constants [6] for Cross Section Formula (10). 


Beam 

Target 

Z 

v, 

y 2 

B 

«1 

V^o 

m 




(mb) 

(mb) 

(mb) 

(mb) 

(GeV 2 ) 

(GeV) 



P 

P 

35.45 

42.53 

33.34 

0.308 

1.0 

5.38 

0.458 

0.545 

P 

n 

35.80 

40.15 

30.00 

0.308 

1.0 

5.38 

0.458 

0.545 


3.1 Proton-Proton cross section 


For proton-proton (pp) cross sections the HZETRN formula works well for pi a b < 1.8 GeV 
and the Particle Data Group parameterizations work well for pi a b >4.7 GeV. Between 
these energies a better fit is obtained with 


a; 


JWN 


7 


pp 


^ 0.16 

Pll 


(13) 
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where 7 = 52.5 mb GeV 0 ' 16 = 158.547 mb MeV 016 . The final new parameterization is 

present 
a pp 

Note that these formulas are discontinuous. See figures 1 and 2. 

3.2 Neutron-Proton cross section 

For proton-proton (pp) cross sections the HZETRN formula works well for pi a b < 1.8 GeV 
and the Particle Data Group parameterizations work well for pi a b >4.7 GeV. Between 
these energies a better fit is obtained with 

<t^ n (T u ) = 40 + 10 cos (a p u - 0.943) exp(-/3 p ° u 8 + 2) (17) 

with pu given by (12). The constants a and /3 are given in Table 3. The final new 
parameterization of the present work is 

present 
a pp 


Note that these formulas are also discontinuous. See figures 3 and 4. 


Table 3: Constants for Cross Section Formula (17). 


a 

(3 

0.00369 MeV” 1 

0.00895741 MeV -0 ’ 8 


°np ZETRN for P 11 < 0.5 GeV (i- e - T u < 124.941 MeV) (18) 

a™ N for 0.5 GeV < p u < 2.0 GeV (19) 

a^ G for p u > 2.0 GeV (i.e. T u > 1271.04 MeV) (20) 


HZETRN 

for 

PP 

JWN 

for 

PP 

n PDG 

a pp 

for 


Pu < 1.8 GeV (i.e. T u < 
1.8 GeV < pu < 4.7 GeV 
pu > 4.7 GeV (i.e. T u > 


1091.74 MeV) 

(14) 


(15) 

3854.69 MeV) 

(16) 
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4 Conclusions 


Nucleon-nucleon total (clastic plus inelastic) cross sections are inputs to space radiation 
transport codes. The nucleon-nucleon total cross sections currently in use in the trans- 
port codes HZETRN and BRYNTRN provide an apparently good global fit to the data. 
However, when one looks more closely at the important GeV energy range the comparison 
to data is poor. This energy range is near the maximum intensity of the cosmic ray spec- 
trum and is therefore important for space radiation studies. Improved parameterizations 
have been developed in the present work, which provide a good fit to data over the entire 
energy range. The new formulas have been incorporated in the HZETRN code. 
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Figure 1: Total proton-proton cross section versus the projectile momentum in the lab 
frame, pi a b- The blue dashed curve is used in HZETRN [4] and the red solid curve is the 
present work. Experimental data are from Reference [6]. 



Figure 2: Detail of Fig. 1. 


13 



Figure 3: Total proton-neutron cross section. The blue dashed curve is used in HZETRN 
[4] and the red solid curve is the present work. Experimental data are from Reference [6] . 



Figure 4: Detail of Fig. 3. 
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